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Abstract  This  paper  reports  an  investigation  of  Computational  Fluid  Dynamics  (CFD)  on  the  influence  of  injec¬ 
tion  momentum  rate  of  premixed  air  and  fuel  on  the  flameless  Moderate  or  Intense  Low  oxygen  Dilution  (MILD) 
combustion  in  a  recuperative  furnace.  Details  of  the  furnace  flow  velocity,  temperature,  CL,  CCE  and  NOv  concen¬ 
trations  are  provided.  Results  obtained  suggest  that  the  flue  gas  recirculation  plays  a  vital  role  in  establishing  the 
premixed  MILD  combustion.  It  is  also  revealed  that  there  is  a  critical  momentum  rate  of  the  fuel-air  mixture  below 
which  MILD  combustion  does  not  occur.  Moreover,  the  momentum  rate  appears  to  have  less  significant  influence 
on  conventional  global  combustion  than  on  MILD  combustion. 
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1  INTRODUCTION 

Increasing  energy  efficiency  has  been  acknowl¬ 
edged  to  be  an  important  issue  related  to  reduction  of 
emissions  of  greenhouse  gases  (e.g.  C02)  and  pollu¬ 
tions  (e.g.,  NOv  and  SOA)  from  industry  combustion  [1, 2], 
Moderate  or  Intense  Low-oxygen  Dilution  (MILD) 
combustion  is  thus  one  of  the  most  important 
achievements  in  combustion  technology  recent  years. 
This  technology,  developed  from  laboratory  tests  to 
industrial  applications,  can  play  a  significant  role  in 
the  mitigation  of  combustion-generated  pollutants 
(NO„  in  particular)  and  greenhouse  gases,  whilst  sat¬ 
isfying  thennal  efficiency  needs.  The  essence  of  the 
technology  is  that  fuel  is  oxidized  in  an  environment 
where  the  oxygen  is  highly  diluted  by  a  substantial 
amount  of  inert  flue  gases  and  temperature  is  above 
that  of  auto-ignition.  Under  these  conditions  tradi¬ 
tional  flames  cannot  be  sustained  and  blow-off  occurs 
due  to  very  strong  shear  motion  caused  by  high  jet 
momentum  and  strong  internal  gas  recirculation. 
Chemical  reactions  occur  in  a  distributed  zone  with  a 
reduced  peak  temperature  [3, 4].  As  a  consequence,  the 
temperature  distribution  is  nearly  uniform,  the  net 
radiation  flux  can  increase  by  as  much  as  30%  [5], 
and  pollutant  emissions,  NOv  in  particular,  are  much 
lower  than  that  from  conventional  flames.  Hence  uni¬ 
formity  of  both  the  temperature  and  the  chemical  spe¬ 
cies  fields  are  distinct  characteristics  of  the  MILD 
technology. 

In  the  literature  this  combustion  is  also  given 
different  names.  Some  relied  on  a  descriptive  form  of 
the  resulting  combustion  process,  i.e.  ‘Flameless  Oxi¬ 
dation’  (FLOX)  or  ‘Colourless  Distributed  Combus¬ 
tion’  (CDC)  since  no  visible  flame  front  develops  [3]. 
If  high  air  preheating  is  used,  the  combustion  is 
named  ‘High  Temperature  Air  Combustion’  (HiTAC 


or  HTAC)  [5].  We  use  the  term  ‘MILD  combustion’ 
throughout  this  paper. 

A  good  volume  of  work  has  been  published  on 
MILD  combustion  and  reviewed  by  Cavaliere  and  de 
Joannon  [4]  and  Tsuji  et  al  [5].  Gaseous,  liquid  and 
solid  fuels  were  investigated  and  some  commercial 
products  have  been  available  for  more  than  one  dec¬ 
ade  [6,  7].  Further  application  of  this  combustion  to 
waste -to-energy  technologies  is  also  underway  [8]. 
Cavaliere  and  de  Joannon  [4]  highlighted  the  current 
knowledge  in  the  field  and  emphasised  on  two  main 
points:  (1)  the  MILD  combustion  regime  is  different 
from  those  of  other  combustions  and  (2)  this  combus¬ 
tion  applies  for  a  narrow  range  of  temperature  that 
allows  design,  optimization  and  adjustment  in  the 
process  by  fine  tuning  external  parameters. 

Although  numerous  fundamental  and  practical 
studies  have  been  perfonned,  fuel-air  premixing  as  a  way 
to  achieve  MILD  combustion  has  received  little  attention. 
Wiinning  &  Wiinning  [3]  noted  that  use  of  premixed 
air-fuel  reactants  under  some  particular  conditions  can 
achieve  this  combustion.  If  the  mixture  can  be  far  off 
stochiometry,  there  will  be  the  opportunity  to  control  the 
temperature  in  a  combustion  chamber  by  the  stoichiomet¬ 
ric  conditions  rather  than  withdrawing  energy  [3]. 

On  the  CFD  modelling  side,  Wiinning  & 
Wiinning  [3]  validated  that  MILD  combustion  can  be 
calculated  with  existing  CFD  simulation  codes.  Christo 
and  Dally  [9]  carried  out  a  detailed  CFD  study  of  jet  in 
hot  co-flow  (JHC)  burner  which  was  designed  to  emu¬ 
late  MILD  combustion  regime.  They  concluded  that  the 
eddy  dissipation  concept  (EDC)  model  with  a  detailed 
kinetic  scheme  offers  best  results.  More  recently,  nu¬ 
merical  modeling  of  MILD  combustion  for  LPG  and 
pulverized  coal  have  also  perfonned  [10-12].  The  flow 
fields  as  well  as  the  temperature  and  oxygen  fields 
were  well  predicted.  Until  now,  MILD  combustion  has 
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been  widely  investigated  by  CFD  and  some  of  numeri¬ 
cal  modelling  results  have  been  verified  by  experiments. 

Our  previous  CFD  work  [13]  has  revealed  that, 
for  non-premixed  fuel  and  air  jets,  there  is  a  critical 
rate  of  the  initial  total  reactant  momentums  below 
which  the  MILD  combustion  does  not  occur.  This  was 
confirmed  by  experiments  in  the  same  study  [13]. 

The  present  study  is  aimed  at  investigating  the 
effect  of  the  momentum  rate  of  inlet  premixed  air  and 
fuel  (natural  gas)  on  MILD  combustion.  The  main 
objective  is  twofold:  (1)  to  identify  the  relationship  of 
flue-gas  recirculation  with  the  momentum  rate  and  (2) 
to  identify  whether  or  not  there  is  also  a  critical  value 
of  the  momentum  rate  present  for  the  premixed  com¬ 
bustion.  The  present  investigation  was  carried  out  us¬ 
ing  the  commercial  computational  software  package 
Fluent  6.3  [14]  as  part  of  our  series  (experimental  and 
numerical)  study  on  the  influence  of  initial  momentum 
rates  of  various  reactant  jets  on  MILD  combustion  in 
different  geometric  configurations  of  furnace  and 
burner  arrangements. 

2  COMPUTATION  DETAILS 

2.1  Furnace  configuration  and  fuel 

The  present  study  simulates  various  premixed 
combustion  cases  occurring  in  a  laboratory-scale  re¬ 
cuperative  furnace,  which  has  been  tested  for  previous 
studies  [13, 15-17].  Here  only  a  brief  description  on  the 
furnace  is  provided  for  the  CFD  domain  establishment 
(Fig.  1).  An  annular  channel  around  a  bluff  body  is 
used  as  a  burner  discharging  premixed  stream  of  fuel 
and  air  into  the  furnace.  The  diameter  of  the  bluff-body 
is  22  mm  while  the  inner  diameter  of  the  tube  varies 
for  generating  different  inlet  momentum  rates  (Table  1). 
Four  exhaust  ports  are  arranged  symmetrically  on  the 
same  wall  with  the  burner  installed.  Two  U-shaped 
cooling  tubes  are  used  to  control  the  heat  load. 

Natural  gas  (90%  CH4  and  10%  C2H6  by  volume 
at  288  K)  is  used  in  this  study.  The  total  volume  flow 
rate  is  about  4x10  3  m3-s  1  and  the  equivalence  ratio 


is  0  =  0.8  .  The  definition  of  equivalence  ratio  is 
0  =  (A/ F)stoic  /(A/ F) ,  where  A  and  F  denote  the  air 
and  fuel  mass  flow  rate,  respectively  [18].  The  air  and 
fuel  used  are  at  ambient  temperature  (288  K). 

2.2  Injection  momentum  rate 

Both  air  and  fuel  are  100%  mixed  together  and 
discharge  as  a  mixture  into  the  furnace  through  the 
annular  tube,  whose  inner  diameter  is  denoted  by  D. 
The  initial  momentum  rate  of  the  mixture  is  defined  as 

Ja  =  |  p0U/dA0 ,  where  p0 ,  Ua  and  Aa  denote  the  mix¬ 
ture  density,  velocity  and  nozzle  exit  area,  respectively. 

Present  calculations  use  six  different  values  of  J0 
(Table  1).  These  values  are  obtained  by  changing  the 
inside  diameter  of  the  burner  D  from  24  mm  to  35  mm. 
As  D  increases,  the  nozzle  exit  area  A0  increases  and 
thus,  to  keep  the  inlet  mass  flow  rate  PQ  =  p0U0A0  con¬ 
stant,  the  injection  velocity  U0  has  to  decrease,  which 
results  in  J0  decreasing.  Table  1  shows  that varies 
from  0.261  kg-m-s  2  (Case  1)  to  0.032  kg-m-s  2  (Case  6). 


Table  1  Case  definition 


Case 

D/mm 

Jet  velocity/nr  s  1 

Jet  momentum/kg-m  s  2 

1 

24 

55.2 

0.261 

2 

26.4 

23.9 

0.113 

3 

28 

16.9 

0.080 

4 

30 

12.2 

0.058 

5 

31 

10.64 

0.050 

6 

35 

6.86 

0.032 

2.3  Computational  conditions  and  models 

The  commercial  computational  software  package 
Fluent  6.3  is  used.  Full  three-dimensional  structured 
grids  are  constructed  with  good  orthogonality.  The 


Figure  1  Schematic  of  furnace  for  CFD  simulation  (in  unit  of  mm) 
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(a)  Computational  domain 
Figure  2 


(b)  Side  view  mesh  (d)  Profile  view  mesh  near 
the  heat  exchanger 

Schematic  diagram  of  CFD  mesh 


final  grids  for  the  half  model  (geometrical  symmetry 
is  utilized)  are  determined  to  be  400000  cells  after 
checking  grid  independence  results  (Fig.  2). 

The  mixture  is  assumed  to  obey  the  ideal  gas  law. 
The  specific  heat  of  the  mixture  is  a  function  of  tem¬ 
perature  (piecewise-polynomial)  [14].  Inlet,  outlet  and 
wall  boundary  conditions  of  the  computational  do¬ 
main  are  set  as  constant  velocity,  constant  pressure 
(static  pressure  =  101  kPa)  and  constant  temperature 
boundaries,  respectively.  The  heat  exchangers  are  de¬ 
fined  as  a  temperature  (400  K)  boundary  and  the  total 
heat  flux  (convection  and  radiation)  through  the  heat 
exchangers  is  checked  after  convergence  (approxi¬ 
mately  3.4  kW). 

The  standard  k-e  model  with  the  standard  wall 
function  is  used  for  modelling  the  turbulent  flow.  The 
turbulent  Schmidt  number  (Sct)  used  is  0.7  [19].  The 
eddy  dissipation  concept  (EDC)  with  two-step  chemi¬ 
cal  kinetic  mechanisms  for  CFI4  and  one-step  mecha¬ 
nism  for  C2FI6  is  applied  [20].  The  present  calculation 
adopts  the  in-situ  adaptive  tabulation  (ISAT)  model  of 
Pope  [21]  to  reduce  the  computational  cost  of  time 
integration.  The  discrete  ordinate  (DO)  radiation 
model  [22]  with  weighted  sum  of  gray  gas  model 
(WSGGM)  to  model  the  gas  emissivity  is  applied  for 
radiation.  The  SIMPLE  algorithm  was  used  for  pres¬ 
sure  velocity  coupling.  The  NO*  formation  is  com¬ 
puted  by  taking  equilibrium  conditions  for  the  reaction 
02^20  [18].  The  formation  rate  is  given  by 

d[NO]/ d?  =  2CNlf  [0]eq[N2] 


kN1{  =  1.8x  1011  exp  (-38370/ J) 

A  second-order  discretization  scheme  is  used  to 
solve  all  governing  equations.  Convergence  is  ob¬ 
tained  when  the  residuals  are  less  than  10  6  for  energy 
and  DO  intensity  and  1 0  1  for  all  other  variables.  The 
outlet  temperature  is  monitored  and  its  variation 
within  1  K  is  allowed  for  convergence  of  the  solution. 

3  RESULTS  AND  DISCUSSION 

Figure  3  shows  contours  of  the  z-component  of 
the  mean  velocity,  v,  (m-s  !),  in  the  central  xz-plane 
(upper)  and  the  xy-plane  at  z  =  400  mm  (lower).  To 
see  the  global  flow  structure  more  clearly,  the  velocity 
vector  field  for  Case  1  is  presented  in  Fig.  4.  Evidently, 
the  furnace  flow  patterns  are  similar  for  all  the  cases: 
i.e.,  they  are  all  characterised  by  a  central  upward  (ini¬ 
tially  fuel-air  mixture)  flow  (v>0),  which  is  reversed 
by  the  furnace  ceiling  (top  wall),  and  side  downward 
flue-gas  flows  mainly  at  four  comers  (v-<0 ).  In  the 
present  furnace  modelled,  the  upward  mass  flow  rate 
( Cup)  must  be  equal  to  the  downward  mass  flow  rate 
(C’down)  at  each  xy  cross-section  above  the  burner  exit, 
i.e.,  Pup(z)  =  Pdowa(z).  So,  the  recirculation  rate  of  the 
flue  gas  can  be  measured  by  the  central  upward  mass 
flow  rate,  viz. 

Pup(z)=  jj  pvz(x,y)dxdy  (1) 

A(z) 

where  A(z)  means  the  area  with  vz>0.  For  the  same 


with 
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Figure  3  Contours  of  the  .--component  mean  velocity  v.  (m  s  ')  in  the  jcr-plane  (y  =  0,  upper)  and  xy-plane  at  -  =  400  mm 

(lower) 


(a)  The  furnace 


(b)  The  central  and  near  field 


Figure  4  Velocity  vectors  in  the  furnace  and  the  central 
and  near  field  (Case  1 :  D  =  24  mm) 


inlet  mass  flow  rate  of  fuel  and  air,  the  higher  value  of 
Pup(z)  reflects  the  more  recirculation  of  the  flue  gas, 
i.e.  PdownCz).  Thus,  Fig.  3  clearly  demonstrates  that  the 
most  dynamical  and  strongest  recirculation  occurs  in 
Case  1,  whose  momentum  rate  Ja  is  highest.  As  J0  de¬ 
creases  from  Case  1  to  Case  6,  Pup  reduces  signifi¬ 
cantly  (and  thus  the  dark-coloured  area  for  contours  of 
vz  ^  6  m-s  1  in  the  cross-section  of  z  =  400  mm 
becomes  smaller).  As  a  result,  the  recirculation  of  flue 
gases  weakens  and,  accordingly,  the  mixing  rate 
between  the  reactants  and  combustion  products 


decreases  greatly. 

The  intense  dilution  of  the  reactants  by  hot  flue 
gas  recirculation  is  considered  as  one  of  the  vital  re¬ 
quirements  for  achieving  MILD  combustion.  The  re¬ 
circulation  controls  the  temperature  and  the  local 
oxygen  concentration  and  hence  the  characteristics  of 
reaction  zone.  Better  mixing  between  hot  products  and 
air  is  required  for  minimizing  hot  spots  and  obtaining 
thermally  unifonn  field,  which  leads  to  lower  levels  of 
NO,  and  CO  emissions. 

The  relative  exhaust  gas  recirculation  rate,  Kv, 
may  be  defined  as  the  ratio  of  mass  flow  rate  of  en¬ 
trained  gases  (i.e.,  the  total  mass  flow  rate  through  the 
effective  jet  cross-section  subtracting  the  mass  flow 
rate  of  inlet  air-fuel  mixture)  to  that  of  the  inlet 
air-fuel  mixture,  i.e.. 


_Pup(z)-Pq(0) 

po(0) 


(2) 


The  upward  mass  flow  rates  (PUp)  and  the  relative  re¬ 
circulation  rates  ( Kv )  for  all  the  cases  were  calculated 
from  their  CFD  data  and  using  Eqs.  (1)  and  (2).  Fig.  5 
shows  the  dependence  of  Kv  on  the  inlet  momentum 
rate  JQ  (kg-m-s  2)  at  five  vy-planes  of  z=  100,  200, 
300,  400  and  500  mm. 

Clearly,  as  the  jet  proceeds  upwards,  i.e.  z  in¬ 
creases,  more  combustion  products  are  entrained  into 
the  upcoming  mixture  jet  and  thus  Kv  is  increased. 
Also  evidently  in  Fig.  5  (a),  increasing  Ja  by  reducing 
D  (nozzle  diameter)  enhances  the  recirculation  of  flue 
gases  and  thus  increases  Kv  at  each  xv-plane.  Such  a 
variation  of  the  .7,, -dependent  recirculation  is  expected 
to  cause  switching  between  the  MILD  and  conventional 


14 


Chin.  J.  Chem.  Eng.,  Vol.  18,  No.  1,  February  2010 


combustion  modes.  Indeed,  this  is  the  case,  as  illus¬ 
trated  by  Fig.  6.  Apparently,  Cases  1-3  have  more 
uniform  temperature  (-1400  K)  distributions  than  do 
cases  4-6  (~2100  K).  It  is  obvious  as  well  that  Cases 
4-6  correspond  to  the  conventional  combustion  mode 
with  visible  flames  in  the  central  region  of  the  furnace 


(a)  At  different  jet  momentum 
z/mm:  ■  100;  O200;  n300;  v400;  «500 


while  Cases  1-3  represent  the  MILD  mode.  Fligh  local 
flame  temperatures  in  the  central  region  of  Cases  4-6 
cause  significantly  more  thermal  NO*  formations, 
relative  to  the  MILD  combustion  of  Cases  1-3  where 
nearly  zero  emissions  of  NO.v  are  produced  (Fig.  7). 
So  it  is  deduced  that  ./„-().()  7  (kg-m-s  2),  the  average 


(b)  At  different  z 

oCase  1;  □  Case  2;  v  Case  3;  •  Case  4;  ■  Case  5;  ACase6 


Figure  5  Exhaust  gas  recirculation  rate  (Kv) 
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temperature/K 


Case  1  Case  2  Case  3  Case  4  Case  5  Case  6 

Figure  6  Temperature  contours  (K)  on  the  central  jrz-plane  (y  =  0) 
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Case  1  Case  2  Case  3  Case  4  Case  5  Case  6 

Figure  7  Contours  of  the  thermal  NO  production  [|rL-(L  dry)"1]  in  the  central  xz-plane  (y  =  0) 
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over  Cases  3  and  4,  can  be  considered  as  the  critical 
value  of  the  inlet  momentum  rate  [indicated  in  Fig.  5 
(a)]  for  this  premixing  arrangement  and  the  furnace 
shown  in  Fig.  1.  Flence,  when  J0<Ja,  the  MILD  com¬ 
bustion  generally  does  not  take  place. 

For  Cases  1-3,  high  recirculation  of  rich-C02 
(~9%)  flue  gas  dilutes  reactants  by  strong  mixing  so 
that  the  volume  (or  mole)  fraction  of  02  is  rapidly  de¬ 
creased  to  a  very  low  level  (~5%,  see  Fig.  8)  before 
reaction  happens.  As  a  result,  when  oxidation  occurs, 
its  reaction  rate  is  considerably  smaller  than  normal 
combustion.  For  Cases  4-6,  low  flue-gas  recirculation 
(ATV  =  0.6— 1.8)  cannot  produce  any  sufficiently  strong 
mixing  of  combustion  products  with  reactants  so  that 
the  initially-high  oxygen  cannot  be  well  diluted  for 
MILD  combustion  to  occur.  On  the  other  hand,  due  to 
low  inlet  speed  and  the  presence  of  a  central  circula¬ 
tion  attached  to  the  bluff  body,  the  upcoming  reactant 
mixture  is  heated  up  so  that  high  temperature  and  rich 
oxygen  enable  oxidation  taking  place  rapidly.  It  is  thus 
seen  in  Figs.  8  &  9  that  traditional  combustion  occurs 


in  the  small  hom-like  zone  attaching  to  the  burner, 
where  all  the  premixed  natural  gas  is  burnt  out  with 
oxygen  to  produce  C02  and  FLO.  Consequently,  the 
volume  fractions  of  02  and  C02  are  nearly  4%  and  9% 
in  the  rest  of  the  furnace  (Figs.  8  &  9).  In  addition,  the 
global  exhaust  CO  emissions  for  Cases  1  to  6  were 
obtained  to  be  20,  20,  2,  5,  21  and  2  pL-(L  dry)  ',  re¬ 
spectively.  This  indicates  that  all  the  combustion 
simulated  is  nearly  complete  no  matter  whether  it  is  in 
the  MILD  or  conventional  regime. 

Furthermore,  several  interesting  observations  can 
be  made  from  Figs.  6,  8  &  9.  Very  small  differences 
occur  between  spatial  distributions  of  temperature,  02 
and  C02  concentrations  for  Cases  4-6.  By  comparison, 
significant  differences  can  be  seen  between  the  same 
distributions  for  Cases  1-3.  It  is  thus  suggested  that 
the  effect  of  the  inlet  momentum  rate  is  less  signifi¬ 
cant  on  conventional  combustion  than  on  MILD  com¬ 
bustion.  Nevertheless,  this  was  not  found  in  our  previous 
investigation  on  the  non-premixed  combustion  [13]. 

Figure  10  shows  where  the  present  Kv  is  located 
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Figure  8  Contours  of  the  mole  fraction  of  02  in  the  central  xz-plane 
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Figure  9  Contours  of  the  mole  fraction  of  C02  in  the  central  jcr-plane 
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amined  by  both  experiments  and  numerical  simula¬ 
tions  [13],  the  premixed  combustion  investigated  pres¬ 
ently  is  also  found  to  have  a  critical  value  of  the  inlet 
momentum  rate  below  which  the  premixed  MILD 
combustion  does  not  occur.  Of  interest  to  note,  as  well, 
the  momentum  rate  imposes  less  significant  influence 
on  the  conventional  combustion  than  on  the  MILD 
combustion.  This  however  may  not  apply  for  the 
non-premixed  combustion  perfonned  in  the  same  fur¬ 
nace,  as  indicated  by  our  previous  investigation  [13]. 
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NOMENCLATURE 


in  the  diagram  of  the  stability  limits  made  by 
Wunning  &  Wiinning  [3]  for  different  non-premixed 
combustion  modes  in  a  furnace  operating  under  ex¬ 
haust  gas  recirculation  mode.  Here  Kv  is  presented 
against  Jav,  the  temperature  averaged  over  the  effective 
jet  cross-section.  This  diagram  quantifies  the  influence 
of  K,  on  the  combustion  process.  For  fK  exceeds  the 
temperature  of  auto-ignition  («800°C),  as  Kv  increases 
beyond  its  critical  value  (depending  on  Tm),  the  com¬ 
bustion  process  develops  from  the  conventional  mode 
to  the  MILD  mode.  This  means  that  Ky  plays  a  sig¬ 
nificant  role  in  determining  the  combustion  mode. 
This  also  indicates  that  the  inlet  air-fuel  momentum 
rate  J0  controls  the  combustion  mode  since  the  mag¬ 
nitude  of  Kv  is  detennined  by  Ja,  see  Fig.  5  (a).  The 
growth  of  J0  enables  Kv  to  increase  beyond  the 
boundaries  of  stable  flame  and  MILD  combustion, 
and  thus  the  combustion  process  switch  from  the  for¬ 
mer  to  the  latter.  However,  the  stability  limit  bounda¬ 
ries  for  the  present  cases  should  be  different  from 
those  obtained  by  Wunning  &  Wunning  [3],  whose 
diagram  may  apply  only  precisely  for  the  non-premixed 
combustion.  Apparently,  most  of  the  data  points  for 
Cases  2-5  are  located  in  the  unstable  flame  region, 
while  those  for  Cases  1  and  6  are  placed  well  within 
MILD  and  stable  flame  regimes  respectively.  Note 
that  the  unstable  flame  cannot  be  simulated  by  the 
present  CFD  work. 

4  CONCLUSIONS 

This  study  has  numerically  investigated  the  in¬ 
fluence  of  the  inlet  momentum  rate  of  the  premixed 
fuel-air  jet  on  combustion  process  in  a  10  kW  recu¬ 
perative  combustor.  The  present  simulations  reveal 
that  the  recirculation  is  strengthened  as  the  inlet  mo¬ 
mentum  rate  increases.  It  is  also  found  that  the  pre¬ 
mixed  MILD  combustion  can  be  established  when  the 
internal  recirculation  of  flue  gases  is  sufficiently 
strong.  Similar  to  the  non-premixed  combustion  ex- 


A  air  mass  flow  rate,  kg-s  1 
A0  the  annulus  tube  exit  area,  m2 
D  outer  diameter  of  inlet  annulus  tube,  mm 
F  fuel  mass  flow  rate,  kg-s  1 

J0  initial  momentum  rate,  kg-m-s  2 

Kv  exhaust  gas  recirculation  rate 

&Nif  reaction  constant,  m3  kmol  1-s  1 
U0  inlet  jet  velocity,  m  s  1 

vz  the  velocity  component  in  the  z-direction,  m-s  1 
p0  density,  kg-m  3 

(j)  equivalence  ratio 
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